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ABSTRACT 

We investigate the nature of tidal effects in compact triple-star systems. The 
hierarchical structure of a triple system produces tidal forcing at high frequencies 
unobtainable in binary systems, allowing for the tidal excitation of high frequency p- 
modes in the stellar components. The tidal forcing exists even for circular, aligned, and 
synchronized systems. We calculate the magnitude and frequencies of three-body tidal 
forcing on the central primary star for circular and coplanar orbits, and we estimate the 
amplitude of the tidally excited oscillation modes. We also calculate the secular orbital 
changes induced by the tidally excited modes, and show that they can cause significant 
orbital decay. During certain phases of stellar evolution, the tidal dissipation may be 
greatly enhanced by resonance locking. We then compare our theory to observations 
of HD 181068, which is a hierarchical triply eclipsing star system in the Kepler field 
of view. The observed oscillation frequencies in HD 181068 can be naturally explained 
by three-body tidal effects. We then compare the observed oscillation amplitudes and 
phases in HD 181068 to our predictions, finding mostly good agreement. Finally, we 
discuss the past and future evolution of compact triple systems like HD 181068. 
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1 INTRODUCTION 

Tidal interactions are known to profoundly impact the orbital evolution of close binary star systems, exoplanetary systems, 
and moon systems. In binary systems, tides drive the components of the systems toward a synchronized state in which the 
orbit is circular, and the components have spins that are synchronized and aligned with the orbit. In a compact triple system, 
no equilibrium state exists, and the endpoint of the orbital evolution is not immediately obvious. Furthermore, the dynamics of 
tidal interactions coupled with multi-body orbital effects can be quite complex and can lead to the formation of astrophysically 
interesting systems (e.g., the Jupiter moon system, short-period exoplanetary systems, and compact binary star systems). 

Although there have been many studies of three-body orbital dynamics including tidal dissipation (e.g., Mazeh & Shaham 
1979, Ford et al. 2000, Eggleton & Kisleva-Eggleton 2001, Fabrycky & Tremaine 2007, Correia et al. 2011), these studies have 
primarily treated tides in a parameterized fashion. Some of these studies take resonant orbital effects into account, but their 
parameterization of tidal interactions ignores resonant tidal effects which may significantly impact the orbital evolution. 
Furthermore, these studies only consider tidal interactions between two components of the system (i.e., they incorporate 
tidal interactions between objects 1 and 2, but ignore interactions with object 3). This approximation is justified in most 
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systems because the third body is relatively distant. However, in sufficiently compact triple systems, a more thorough study 
of three-body tidal effects is necessary. 

Observing tidal effects is very difficult because of the long time scales associated with tidal evolution, and tidal orbital 
decay has only been observed in a few rare circumstances (such as the orbital decay of a pulsar-MS binary, see Lai 1996, f 997; 
Kumar & Quataert 1998). Until recently, the direct observation of tidally induced oscillations was difficult because of the 
extreme precision required (~ one part in a thousand) over base lines of several days. Fortunately, the continuous observation 
and high accuracy of the Kepler satellite is allowing for the direct observation of tidal effects (e.g., Welsh et al. 2011, 
Thompson et al. 2012) and detailed analyses of tidally excited stellar oscillations (Fuller & Lai 2012, Burkart et al. 2012). 
Recently, observations of luminosity variations in the compact triple system HD 181068 (also known as KIC 5952403 or the 
Trinity system, see Derekas et al. 2011) have provided evidence for three-body tidal effects. To our knowledge, no detailed 
study of dynamical tides exists for compact triple systems. 

We present the first detailed investigation of the tidal excitation of stellar oscillation modes in stars in compact triple 
systems. Hierarchical three-body systems can experience tidal forcing at frequencies unattainable for two-body systems, 
allowing for the tidal excitation of high frequency p-modes in the convective envelopes of the stellar components. We investigate 
the observational signatures of three-body tides, calculate the amplitudes to which stellar oscillation modes are tidally excited, 
and study the orbital evolution induced by three-body tides. 

We also compare our theory to the observations of luminosity fluctuations in HD 181068, accurately characterized by 
Derekas et al. 2011 and Borkovits et al. 2012. HD 181068 is a triply eclipsing hierarchical star system, with a red giant 
primary orbited by two dwarf stars. The dwarf stars orbit each other every 0.906 days, and their center of mass orbits the 
primary every 45.47 days. In addition to eclipses, the lightcurve shows oscillations at combinations of the orbital frequencies. 
We demonstrate that these oscillations are tidally excited oscillations in the red giant primary excited by the orbital motion 
of the dwarf stars. 

Our paper is organized as follows. In Section [2] we derive the strength and frequencies of tidal forcing unique to three- 
body systems. In Section[3]we calculate the amplitudes to which modes are tidally excited, and estimate the resulting stellar 
luminosity variations. In Section [4] we calculate the orbital evolution induced by the tidally excited modes. In Section [5] we 
describe observations of tidally excited modes in HD 181068, and in Section [6] we compare these observations to our theory. 
In Section [7] we calculate the possible past and future orbital evolution of systems like HD 181068. Finally, in Section [8] we 
discuss our results. 



2 THREE-BODY TIDAL FORCING 

Let us consider a triple-star system composed of a central primary star (Star 1), orbited by a pair of companion stars (Stars 
2 and 3) at frequency Qi. Stars 2 and 3 orbit one another at a much higher frequency fl23 3> fii. Figure [I] shows a diagram 
of the orbital configuration. We adopt a coordinate system centered on Star 1 with z— axis in the direction of its spin angular 
momentum vector. We consider the case in which all three stars have circular coplanar orbits that are aligned with the stellar 
spins. We choose a direction of reference such that the observer is located at the azimuthal angle <f> = 0. 

We wish to calculate the form of the time-varying gravitational potential of the short-period binary (Stars 2 and 3) as 
seen by Star 1. The tidal potential due to Stars 2 and 3 can be decomposed into spherical harmonics as 

4tt 



Ulm = 
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Here, M2, D2, 62, and cf>2 are the mass, distance, polar angle, and azimuthal angle of Star 2 (M3, D3, #3, and <j)3 are the same 
quantities for Star 3), while G is the gravitational constant. The dominant terms have I — \m\ = 2 and I = 2, m = 0, so we 
will consider only these terms in our analysis. Since we restrict our analysis to coplanar orbits aligned with the primary spin, 
#2 = 3 = it/2. 

We wish to express the potential in terms of the angular orbital frequencies f2i and fi23- Some trigonometry reveals that 
the distance from Star 1 to Star 2 is 

D % = a * + (m^Tm^ 23 ) 2 + 2 M2TM3 aia23COsW23 " (2) 

where (j>i is the phase of the orbit of Stars 2 and 3 about Star 1, $23 is the phase of the orbit of Star 2 around Star 3, ai is 
the semi-major axis between Star 1 and the center of mass of Stars 2 and 3, and a23 is the semi-major axis between Stars 2 
and 3. Also, the law of sines reveals 

sin(<fe - <t>i) = M2 M + 3 Ma g sin(fe - (3) 
In a hierarchical triple system, (123 < ai, so we expand distances and angles in powers of the small parameter 

E5023/01. (4) 
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Figure 1. This diagram (not to scale) depicts the geometry of a triple-star system at arbitrary orbital phase. a\ is the semimajor axis 
between the center of Star 1 and the center of mass of Stars 2 and 3. 0,23 is the semimajor axis between Stars 2 and 3. D2 is the distance 
between Stars 1 and 2. <f>\ is the orbital phase (relative to a direction of reference) of the center of mass of Stars 2 and 3, while (j>2 is the 
orbital phase of Star 2 about Star 1. if>23 is the orbital phase of Stars 2 and 3 about each other relative to the same direction of reference. 
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The values of D 3 3 and e~ lm< ^3 are the same as equations ^ and (JsJ) , but with 2 and 3 subscripts reversed, and with 

023 — > 023 + 7T. 

Inserting equations ([5| and |6| and their counterparts for Star 3 into equation we find that the zeroth order contri- 
bution (the term proportional to <r) for I — 2 to the tidal potential is 



r/(0) _ 
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where W2,o = —\pxi§ and W2,±2 = \/37r/10. To zeroth order, the tidal potential due to stars 2 and 3 is simply that of a 
single star of mass M2 + M 3 at a semi-major axis ai. 

Upon summing the contributions from both stars, the first-order terms (proportional to e) in equation ([T]) vanish. However, 
the second-order terms produce an additional component of the tidal potential. Dropping terms with no 0i or 023 dependence 
because they have no time dependence, we find that the second-order component of the tidal potential is 
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Then the m = component yields 
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The first term in brackets produces retrograde forcing, i.e., it will excite modes with negative angular momentum that rotate 
in the opposite direction to the orbital motion of Stars 2 and 3 about one another. The second term in brackets produces 
prograde forcing, exciting modes with positive angular momentum. 

We define 0i — Qit and 023 = 02, o + f223t, where 02, o is the phase 2 3 of the short-period orbit t — 0, when their center of 
mass is closest to the line of sight. The second-order m — component of the tidal potential produces forcing at a frequency 
(To = 2(f2 23 — fii). The m — ±2 components produce forcing at two frequencies, cr_ 2 = 2f2 2 3 and 02 = 2(0.23 — 2f2i). These 
frequencies arise because they are the conjunction frequency in a frame corotating with the long-period orbit, but are Doppler 
shifted by mfii when viewed in the inertial frame. Due to the hierarchical structure of the system, these frequencies can be 
larger than the dynamical frequency of Star 1, causing tidal forcing at high frequencies seldom occurring in binary systems. 

Equation (JsJ) describes the tidal potential in the inertial frame, but we also wish to calculate the forcing frequencies in 
the rotating reference frame of Star 1. To do so, we make the transformation = 0, + f2 s i, where r is the azimuthal angle 
in the rotating frame and Q s is the angular spin frequency of Star 1. Then in the rotating frame, 
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Equations ( 13 1— ( 15 1 are general for circular coplanar orbits aligned with the stellar spin. If Star 1 is synchronized with 
the orbit of Stars 2 and 3 such that Q s> i — fiiQ 
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Thus, when the spin of Star 1 is aligned with the orbit of Stars 2 and 3, the absolute values of the forcing frequencies of all 
three modes (the axisymmetric m = mode, and the prograde and retrograde m = ±2 azimuthal modes) in the corotating 
frame are identical, namely, \u\ = 2(f223 — fii). 



3 MODE EXCITATION AND OBSERVATION 
3.1 Mode Amplitudes 

With the tidal potential known, we may calculate the amplitude of the tidally forced modes. In this section we consider only 
tidal forcing by the I = 2 second-order component of the tidal potential. We use the subscript a to refer to a mode of angular 
degree I and m, with a frequency in the rotating frame of uj a . In this frame, the mode amplitude c a satisfies (Schenk et al. 
2001) 

<i 

c a + (iuj a +j a )c a = - — (£ a (r),-Vf7) 

_ iAgQg i m (n s -f2i)t 
2e a 



F„ 



O 2i(fi 2 3-«l)* + 2i<£2,0 



+ F- 



-2s(n 2 3-fil)t-2i^2,0 



(17) 



1 In compact triple systems in which three-body tides operate, it is very likely that f2 s i ~ fli because the tidal synchronization time 
due to the zeroth order tidal potential of equation [7] will be much shorter than tidal time scales due to the three-body tides of equation 
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Here, ui a = e a — mC a ^sSgn(e a ) is the mode frequency in the rotating frame, e a is the unperturbed frequency for a non- 
rotating star, C a is the rotational kernel (see Fuller & Lai 2012), and 7 a is the mode damping rate. In this formalism, both 
m and e a can have positive or negative values. The dimensionless coefficient Q a is the overlap integral of the mode with the 
tidal potential, defined as 



where £ a is the mode eigenvector normalized via the condition 

(UU = m 1 r[, 

and Ri is the radius of Star 1. 

The non-homogeneous solution of equation ( 17 1 is 

A a Qa 
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where 

I / m = 2(n 2 3-n 1 ) + m(n s -ni). (21) 

The total tidal response of the star is £(r,t) — c a (i)£ a (r). We assume that the star is slowly rotating (flj, <C uj a ) such 
that £ a (r) and Q a are independent of m or Q, s , and u) a ~ e a . Summing over both signs of uj a yields 
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Summing over both signs of m then yields 
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The m — modes produce a radial displacement 

£ r (r, t, m = 0) = ^ 2AaQaF m Da,m£r,a(r)P 2 fi(cOs9) COs(u m t + 202,0 + t/j a ,m)- 
a ,oj >0,m — 

The m = ±2 modes produce a radial displacement 
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The first term is the radial displacement due to the retrograde mode, and the second is due to the prograde mode. 



(24) 
(25) 
(26) 
(27) 



(28) 



3.2 Luminosity Variations 

The luminosity variation produced by a mode, AL a /L, is not trivial to calculate, because the temperature fluctuation produced 
by the mode is sensitive to non-adiabatic effects and subtleties in the outer boundary condition. We refer the reader to other 
studies (e.g., Buta & Smith 1979, Robinson et al. 1982) which have attempted to quantify the visibility of a mode given its 
amplitude and eigenfunction. If the mode is adiabatic, its luminosity fluctuation is 

—j^- = 2AaQaFmDa, m £r,a(R)VaYl im (Q , O ) COs((7 m t + 202,0 + ^a,m) (29) 

with 



a m = 2(ri 2 3 — fli) — mfii, 
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and here m = 2 corresponds to the retrograde mode and m = —2 corresponds to the prograde mode. The angular coordinates 
9 and <j> indicate the direction of the observer relative to the symmetry axis of the mode, here assumed to be the spin axis of 
Star 1. We define V a to be the visibility function of the mode, which is dependent on geometrical, opacity, and non-adiabatic 
effects. In high-inclination systems like HD 181068, 9 ~ 90°, and we have chosen our coordinate system such that <j>o = 0. In 
this case, Yi >m (0 o , (f> a ) ~ — \/5/(167r) for m = modes and Yi, m (# , 4>o) — \/l5/(327r) for m = ±2 modes. 

The value of V a is uncertain. If non-adiabatic effects are significant, V a will be complex and the luminosity variation will 



be phase-shifted from equation (29 1. Here we consider adiabatic modes, but we remember that an observation of a phase shift 



is indicative of non-adiabatic effects. In the adiabatic limit, Buta & Smith (1979) find 



ai + ft + 4V ad 7,i/ Q 



(31) 



where ai, /3i, and 7; are coefficients of order unity that depend on the spherical harmonic I of the mode and the stellar 
limb-darkening function (we use the values given in Buta & Smith 1979), and V a( j is the adiabatic temperature gradient at 



the surface. The first term in equation (31 1 is due to surface area distortions, the second is due to surface normal distortions, 



and the third term is due to temperature effects. The function H a describes the magnitude of the pressure perturbation 
compared to the radial surface displacement, i.e., AT a /T = H a AP a / P. According to Dziembowski (1971), the value of H a 
for adiabatic oscillations is 



Ha = 



1(1 + 1) 



-4- 



where 



Wdyn 

In this description, H Q 



(32) 



(33) 



is large for high-order g-modes (low 



and ojdyn = \j GM\ / Rf is the dynamical frequency of Star 1. 
values of ui a ) or high-order p-modes (high values of ui a )- 

Gouttebroze & Toutain (1994) have attempted to estimate the value of V a for p-modes modes in a solar model. They 
calculated non-adiabatic mode eigenfunctions, and then calculated luminosity variations by adding up the perturbed flux from 
a grid of emitting surface elements. They found that V a is of order unity for low-order p-modes, but has substantially larger 
values for the high-order p-modes that typically produce solar-like oscillations. In the analysis of Section [6j we calculate V a 
from equation ( |31[ ), considering a lower limit of H a ~ 1 and an upper limit of equation ( |32[ ). 



4 EFFECT OF MODES ON ORBITAL EVOLUTION 
4.1 Hamiltonian Formalism 

The tidally excited modes draw energy and angular momentum from the orbital motions and deposit them in Star 1. The 
effect of the modes on the orbits can be calculated from the Hamiltonian of the gravitational interaction between the modes 
and the stars. This Hamiltonian is 



H = j d\U{r,t)Y,c a {t)Spl{v) 



(34) 



where 8p is the Eulerian density perturbation associated with each mode. Performing the integration over the volume of Star 
1, and considering only / = 2 terms yields 
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Inserting equation (201 into equation (|35h , we find 



H = -MrRl [A a Q a ) 

a,u>>0 



2^Q 



F m e 



+ 



2i(023-^l)-* m ^l I T? _-2j(023-</'l)- im </ > l 

i 1 — mt; 



(36) 



with a fi — 2(^23 — Jli) — Wifli, and a/2 = —2(^23 — fii) — mfii. 

The orbital evolution equations are L\ = —dH/d(f>i, L23 = —dH/d(f)23, where L\ = fiiaffli and L23 = ^23(123^23 are the 
angular momenta of the outer and inner orbit, respectively, and /xi — M\(Mi + M3)/ (Mi + M2 + M3). In the limit Q a <C uj a 
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we find 
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We have discarded rapidly oscillating terms because they produce no secular variations. The torque due to the m = modes 
is 



L 1 (m = 0) = 8M 1 Bi ^ {A a Q a ) 2 FlY a , 



where 



r — 
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and L2i(m = 0) = —Li(m = 0). The torque due to m = ±2 modes is 
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(40) 
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Equations ( 39|42 I may be written in terms of a dimensionless damping rate T a , m — r Q , m w| yn , such that we obtain the 
more familiar scaling for tidal dissipation: 



it n\ o GMi 2 / A*23 
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for m = ±2 modes. 



Conservation of angular momentum requires that the torque on Star 1 is L„ = —L\ — L23, or 



L, = 



GM 2 2 /^23 

Ri ' Ui 
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We can see that m — modes draw angular momentum from the inner orbit and transfer it to the outer orbit. The m — 2 
(retrograde) modes take angular momentum from the inner orbit and the spin of Star 1, depositing it in the outer orbit. 
The m = —2 (prograde) modes draw angular momentum from the inner orbit and transfer it to the spin of Star 1. Thus, in 
all cases, three-body tides cause the inner orbit to decay and the outer orbit to expand, although the outer orbit expands 
by a much smaller factor because of its larger moment of inertia. Star 1 can be either spun up or spun down, depending 
on which mode (prograde or retrograde) contains more energy. Under most circumstances, Star 1 will usually be spun down 
because F 2 > F- 2 . In the limit of zero eccentricity, the orbital energies change as E\ = SI1L1 and -B23 = £123^23 ■ Consequently, 
ei = e23 = 0, and the orbits remain circular. Finally, we note that three-body tides tend to increase the stability of hierarchical 
triples because they cause the period ratio P1/P23 to increase. 
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Table 1. Properties of the HD 181068 system as measured by Borkovits et al. 2012. 





M(M ) 


R(Rq) 


T eS (K) 


Star 1 


3.0 ±0.1 


12.46 ±0.15 


5100 ± 100 


Star 2 


0.915 ± 0.034 


0.865 ±0.010 


5100 ± 100 


Star 3 


0.870 ± 0.043 


0.800 ±0.020 


4675 ± 100 



Table 2. The significant peaks of the period analysis for HD 181068. A number of peaks located < O.ld 1 wore left out of the analysis. 



No. 


Frequency 


Orbital Relation 


Amplitude 


Phase 


S/N 




(d- 1 ) 




(mmag) 






h 


2.1203 


2(Q 23 -2U 1 )/(2tt) 


0.44 


0.0403 


58 


h 


2.1643 


2(fi 23 - fil)/(27T) 


0.25 


0.2889 


35 


h 


1.1065 


n 23 /(2TT) 


0.08 


0.2843 


3.4 


h 


2.2084 


2n 23 /(27r) 


0.048 


0.8719 


6.8 



5 OBSERVATIONS OF HD 181068 

We now apply our theory to the compact triple-star system HD 181068, also known as the Trinity system, whose properties 
are listed in Table 1 (see Derekas et al. 2011, Borkovits et al. 2012). Stars 2 and 3 orbit about Star 1 at an angular frequency 
f2i = 0.138 d -1 , while Stars 2 and 3 orbit about each other at fl 23 = 6.94 d _1 . All three stars are nearly exactly coplanar, 
with orbital inclinations of i ~ 90°. The red giant primary has a radius and surface temperature typical of red clump stars, 
with Rt = 12.46 ± 0.15i?© and T cfi = 5100 ± 100K, but has a fairly large mass of Mi = 3.0 ± O.1M . 

In the discovery paper of HD181068 (Derekas et al. 2011), 218 days of Kepler data were analyzed. One of the most 
surprising results was that the main component of the system, which is a red giant star, did not show solar-like oscillations. 
Instead, other pulsations were detected with frequencies close to double the orbital frequency of the short-period binary. 
Derekas et al. 2011 suggested that these pulsations might be tidally forced oscillations and that there might be a mechanism 
which suppresses the solar-like oscillations. 

5.1 Fourier Analysis 

As of writing, 11 quarters of data from Kepler have been made available to us, representing almost three years of essentially 
uninterrupted observations. The first six quarters of data were obtained in long-cadence mode (one point every 29.4 minutes), 
while Q7 to Qll were taken in short-cadence mode (one point every 58.9 seconds). Here we present a period analysis of the 
combined long-cadence and short-cadence data, the latter covering 450 days. 

As a first step, we subtracted the eclipses and rotational variations by using the light curve fit of Borkovits et al. 2012, 
which resulted in a nearly continuous data set containing the pulsations. For the period analysis, we used Period04 by Lenz & 
Breger 2005. In the Fourier-spectra, a number of the peaks were located < O.ld -1 indicating long term variability, remnants 
of the light curve fit, or instrumental effects that are negligible in the present analysis. Peaks with significance are listed in 
Tableland the Fourier-spectrum is shown in Figure [2] 

The most intriguing result of the period analysis is that /i and / 2 are linear combinations of the two orbital frequencies, 
suggesting a tidal origin. Their angular frequencies are separated by 20i , creating a beat pattern in which they are in phase 
near the primary eclipses and occultations (see Figurej3|. The peaks frequencies / 3 and correspond to one and two times the 
orbital frequency of Stars 2 and 3, and may be caused by the imperfect subtraction of the eclipses or spots on the components. 
However, they could have a tidal component as well. 

5.2 Lack of solar-like oscillations 

Red giant stars are well known to show solar-like oscillations (De Ridder et al. 2009, Chaplin et al. 2011a) which are excited 
by near-surface convection (see, e.g., Houdek et al. 1999, Samadi et al. 2007). Nearly all red giants observed by Kepler show 
detectable oscillations (Huber et al. 1010, Kallinger et al. 2010, Hekker et al. 2011), which can be used to probe the internal 
properties of the star (e.g., Bedding et al. 2011, Beck et al. 2012, Mosser et al. 2012). The frequency of maximum power 
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Figure 2. Fourier spectrum of 11 quarters of long-cadence data. 
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Figure 3. Red curve: Sample of light curve of HD 181068 between long duration eclipses. The x— axis is the phase of the long-period 
orbit, measured from the primary minimum at BJD 55545.466. Green curve: The simulated light curve of the oscillations. Note the 
beating pattern due to two close frequency oscillations. 



(I'max) of solar-like oscillations can be estimated using the scaling relation (Brown et al. 1991) 

_ M/M e (T cB /T cB , e f- 5 

fmai — ^77 fmax.Q, {lO ) 

L/Lq 

with t'max.e = 3090/^tlz. Using the fundamental properties given by Borkovits et al. 2012, we calculate f ma x = 64 ± 7/iHz 
(5.5 ± 0.6 c/d) for the red giant component. To search for solar-like oscillations in HD181068, we first remove all primary 
and secondary eclipses from the data using the ephemeris given by Borkovits et al. 2012. To further remove long-periodic 
variability we then apply a Savitzky-Golay filter (Savitzky & Golay 1964) with a width of 5 days. Finally, we pre-whiten the 
most significant low frequency variations due to tidal oscillations, as discussed in the previous section. 

The upper panel of Figure [4] shows a power spectrum of the residual Ql-Qll long-cadence light curve. For comparison the 
middle panel shows the power spectrum of KIC4662939, a Helium-core burning red giant with similar fundamental properties 
as red giant component in HD181068 (Bedding et al. 2011). Despite an increase of the data set length by a factor of three 
compared to Derekas et al. 2011, the lack of solar-like oscillations in HD181068A is clearly confirmed. On the other hand, 
we observe that both stars exhibit a similar decrease of power from low to high frequencies, which is the typical signature 
of granulation (see Mathur et al. 2011). This confirms that both stars indeed have similar fundamental properties, but that 
solar-like oscillations are suppressed in HD181068A. 

We speculate that the lack of solar-like oscillations of the main component is related to the close multiplicity of the 
components of HD 181068. Derekas et al. 2011 measured a rotational velocity of Star 1 of v sini = 14±1 km s _1 , corresponding 
to Q a i = 0.14±0.01 d" 1 . This spin frequency is consistent with Star 1 being synchronized with the long-period orbit, indicating 
that it has been tidally synchronized. It is also an abnormally large spin frequency for red giants (de Medeiros et al. 1996), 
and the rapid rotation may generate a strong magnetic dynamo. Indeed, the light curve of HD 181068 exhibits flaring events 
(some events occur during occultations of Stars 2 and 3, indicating the flares originate from Star 1, see Borkovits et al. 2012) 
that indicate high levels of magnetic activity. Chaplin et al. 2011b has shown that solar-like oscillations are suppressed in 
abnormally active stars, presumably due to their rapid rotation or due to the effects of strong magnetic fields. We therefore 
speculate that the tidal synchronization of the primary in HD 181068 creates rapid rotation and high magnetic activity that 
suppress the excitation of solar-like oscillations. 

The available short-cadence data also allow us to search for solar-like oscillations in the dwarf components. Using equation 
(48 1, we expect the dwarf components to oscillate at frequencies between 3500 — 5000/iHz. The bottom panel of Figure [4] 
shows the power spectrum of the short-cadence data after performing the same corrections as described above. We do not 
detect any significant power excess in the data. Given the much lower luminosities of the dwarf companions and the amplitude 
dilution by the brighter giant component, this implies that the amplitudes are too small to be detected with the data at hand. 
Solar-like oscillations in the dwarfs may be also suppressed by the same mechanisms described above. 
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Figure 4. Top panel: Power spectrum of the Ql-Qll long-cadence data of HD 181068 after removing eclipses and low frequency 
variability. The red line shows the power spectrum smoothed with a Gaussian filter with a full-width at half maximum of 5/xHz. Middle 
panel: Same as top panel but for KIC 4662939. Bottom panel: Power spectrum of the Q7-11 short-cadence data of HD 181068 after 
removing eclipses and low frequency variability. 



6 COMPARISON WITH OBSERVATIONS 
6.1 Stellar Model 

We generate a M = 3.OM , R = 12.4R , T cS = 5100K, z = 0.015 helium burning red giant stellar model using the MESA 
stellar evolution code (Paxton et al. 2011). Figure[5]shows a propagation diagram for the stellar model. The high Brunt- Vaisala 
and Lamb frequencies in the radiative interior of the star allow high-order g-modes to propagate in the stellar interior, while 
the small Brunt- Vaisala and Lamb frequencies in the convective envelope of the star allow for p-mode propagation. 

Using the model described above, we calculate the adiabatic stellar oscillation modes using the usual boundary conditions 
(see e.g., Unno et al. 1989), but normalized according to equation | |19[ ). Figure [6] shows the values of £, r ,a{R) and Q a as a 
function of u) a for our stellar model. Low frequency modes g-modes (lower than uj a « 1) are trapped in the core. Higher 
frequency modes still have g-mode character in the core, but have p-mode character in the convective envelope of the star. 
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Figure 5. Top: the Brunt Vaisala frequency N (solid line) and Lamb frequency L2 (dashed line) as function of radius r. The red 
horizontal (overlapping) lines mark the eigenfrquencies Q a of the modes shown in the bottom panel. Bottom: mode eigcnfunctions £ r , a 
for an envelope mode with u a = 6.66 (long dashed line), a neighboring mixed mode with Q a = 6.71 (dotted line), and a core mode with 
Q a = 7.17. All quantities are calculated for our M = 3.0Mq, R = 12.4i?0 helium burning red giant model, and are plotted in units of 
G = M = R = 1 . The modes are normalized via equation ||19| . 



Modes whose inertia lies primarily in the convective envelope are the envelope dominated mixed modes, while the inertia of 
the neighboring mixed modes is split between the core and envelope. Figures [5] and [6] demonstrate that the envelope modes 
have large values of £ r ,a(R) and Q a compared to neighboring mixed modes. Consequently, the envelope modes are easily 
excited to large amplitudes (because of the larger values of Q a ) and produce large luminosity fluctuations [because of the 
larger values of £ r ,a(R)]- The envelope modes thus dominate the visible and energetic response of the star to high frequency 
tidal forcing. 



6.2 Comparison with Observed Luminosity Fluctuations 

Here, we compare our theory to observations of HD 181068. Unfortunately, it is difficult to predict the amplitude of a tidally 
excited mode, even if the system parameters are relatively well constrained. Part of the reason is that the value of the frequency 
detuning -D Q , m (see equation 251 can vary by orders of magnitude over small changes in uj a . Since the values of ui a depend 
on the precise mass, radius, and internal structure of the star, accurately calculating D a>ra for each mode is very difficult. 
Instead, we choose to calculate the amplitude of the luminosity fluctuations as a function of the stellar radius Ri , which is 
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Figure 6. Radial displacement at the stellar surface £ r ,a(R) (top) and overlap integral Q a (bottom) for stellar oscillation modes as a 
function of mode eigenfrequency uj a . The modes are calculated from our M = 3.0Mq, R = 12 ARq red giant model. The modes at the 
tips of the spikes are envelope dominated mixed modes, while the surrounding modes are mixed modes whose inertia is split between the 
core and envelope. The vertical lines mark the values of 2f2i and 2f223 in the HD 181068 system. 



constrained to be R\ = 12.46 ± O.15i?0. We compute the dimensionless values of 0J a for the stellar model described in Section 



6.1 and then scale them to dimensional frequencies uj a using equation (33 1. The range of radii R\ is meant to encompass 



uncertainties in u> a , Mi, R±, 0. s ,i, non-adiabatic effects, etc., that affect the precise values of ui a in HD 181068. 

Figure[7]displays the theoretical and observed luminosity variations produced at the tidal forcing frequencies 2(^23 — 2f2i), 
2(^23 — fii), and 2S723- Recall that the observed variation at 2O23 may be contaminated by imperfect eclipse subtraction 
and/or spotting effects. The luminosity variations are calculated using equation ( 29 1, using Q23 and fli observed in HD 181068, 
and using M2 = M3 = Q.QMq and o — 87.5° (as measured by Borkovits et al. 2012). The top panel uses H a = 1, while the 
bottom panel uses The luminosity variations have peaks and dips at values of Ri for which the 

frequency of a mode is nearly resonant with the forcing frequency. For H a — 1, the theoretical luminosity variations are well 
below the observed variations, except very near a resonance with an envelope p-mode (i.e., at R 1 ~ 12.8R@). Unless the HD 
181068 is in a resonance locking state (see Section [?b, it is unlikely to observe it so close to resonance. 

In contrast, the theoretical luminosity variations are closer to the observed variations for H a = [1(1 + 1)/oj 2 ~ 4 — ui%] . The 
best agreement is obtained for Ri ~ 12. 6Rq, and matches the luminosity variation for each theoretical and observed frequency 
to within a factor of 2. We conclude that the luminosity variations in HD 181068 are likely dominated by temperature effects, 
characterized by large values of H a . The ordering of the amplitudes of the oscillations is naturally explained by our theory, 
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Figure 7. Luminosity variations AL/L as a function of stellar radius, using H a = 1 (top panel) and H a = u) a (bottom panel), for our 
stellar model. The actual stellar model has R\ = 12.4i?0, but to make this plot we consider identical models scaled to different radii. 
The black lines are the luminosity variations of the mode at frequencies 2(f!23 — 2fii) (dashed line), 2(Q23 — Sli) (solid line), and 2Q.2Z 
(dotted line). The horizontal red lines are the observed luminosity variations in HD 181068 at the same frequencies. 



i.e., AL/L(m — 2) > AL/L(m — 0) > AL/L(m = —2), because F% > F > F_2. However, the measured amplitude ratios are 
slightly different than what we expect away from exact resonances, and the cause of the discrepancy is unclear. 

For non-adiabatic modes, our theory also predicts the phases of the observed oscillations. In particular, equation (29 1 
shows that for non-resonant adiabatic modes (i.e., V a is real and ip a — or ij) a — it), the phase difference between two 
oscillation frequencies is A<f> ~ mcj>i. This implies that the prominent frequencies /1 and /2 (which have m = 2 and m = 0, 
respectively) should be in phase at fa ~ and fa ~ 7r, i.e., they should be in phase during the eclipses and occultations of 
Star 1. Indeed, the measured oscillations are in phase at these times, as can be seen from the simulated light curve in Figure 
[3] This suggests that the observed luminosity fluctuations are being produced by non-resonant modes, consistent with our 
findings above. 
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7 ORBITAL EVOLUTION OF HD 181068 

We now wish to calculate the orbital evolution induced by the tidally excited modes in realistic systems such as HD 181068. 



The quantity in the orbital evolution equations ( 43 1-( 46 1 with the most uncertainty is the damping rate j a . An accurate 
calculation of 7 a requires a fully non-adiabatic calculation of mode eigenfrequencies, and must also include the turbulent 
damping of modes in the convection zone, which is not well understood. We choose to estimate y a using a quasiadiabatic 
WKB damping rate (see Burkart et al. 2012). We check that the damping rate calculated in this manner is the same order 
of magnitude as the damping rate inferred from observations of mode lifetimes of solar-like oscillations in red giants (see 
Belkacem 2012) for modes with uj a ~ v u 



'max ■ 



To understand the effects of tidally excited modes on the orbital evolution of a realistic triple system, we calculate 



the orbital evolution of a system resembling HD 181068 using equations ( 43 1-( 46 1. We consider coplanar, circular orbits, as 
observed in HD 181068. The orbital frequencies change as 

fii=-^ (49) 

and 

n 23 = -^r, (50) 

^23^3 

and the orbital semi-major axes change as a/ a — — 2f2/(3f2). We define the tidal dissipation time scale as 

_ a 2 3 _ M23«23^23 , ci s 

Itidc — - — — ■ ■ (.0 rj 

«23 21/23 

The orbital evolution timescales due to tidally excited modes can be comparable to the stellar evolution time scales of 
the stars in the system. Hence, it is important to evolve the properties of the stars in the system simultaneously with the 
orbital elements. In our evolutions, we compute the radius of Star 1 as a function of time using the MESA stellar evolution 



code (Paxton 2011). The changing value of Ri affects not only the values of L in equations ( 43 (- ( 46 ( , but it also affects 



the eigenfrequencies u) a because the stellar oscillation frequencies scale as ui a oc GMi / Rf . Accounting for the changing 
eigenfrequencies is important because it can lead to resonance locking (Witte & Savonije 1999, Fuller & Lai 2012), allowing 
for tidal evolution on stellar evolutionary timescales rather than the longer non-resonant tidal evolution time scales. 

In our calculations, we do not calculate new eigenfreqencies Q a and associated eigenfunctions at each time step. We find 
that in the red giant phase of stellar evolution, the variations in Q a as a function of time are relatively small, and that the 
variations in u: a are dominated by variations in ojdyn due to the changing value of -Ri . We use the values of Q a and Q a shown 
in Figure [6] and calculate ui a from equation ( 33 l . 



In addition to the tidally excited modes discussed in this paper, our evolutionary calculations should account for "two- 
body" tidal effects, i.e., the tidal effects due to the zeroth order component of the tidal potential in equation Q. These tidal 
forces have no dependence on the small parameter e, and hence they act on much shorter time scales. Furthermore, tidal forces 
between Stars 2 and 3 will act on even shorter time scales due to their close separation. Therefore, two-body tidal forces cause 
Stars 2 and 3 to have a circular orbit around each other, rotating synchronously with that orbit. They also cause the center 
of mass of Stars 2 and 3 to have a circular orbit around Star 1, with Star 1 rotating synchronously with that orbit. Thus, in 
our evolution, we enforce = Q,\ and f2 s ,2 = ^s,3 = at all times. We account for the angular momentum redistribution 
associated with these processes, although the stellar spins contain only a small fraction of the total angular momentum of the 
system. 

Because the moment of inertia of Star 1 is much less than that of the orbit of Stars 2 and 3 about Star 1, a small amount 
of orbital angular momentum deposited in Star 1 by three-body tidal effects can drastically change its spin frequency. This 
angular momentum will then be transferred back to the orbit of Stars 2 and 3 about Star 1 by two-body tidal effects until 
synchronism is restored. Therefore, after each time step in our evolution, we calculate ALi, AL23 and AL* from equations 



( 43 (- ( 47 1 . We then adjust the values of Q. s ,i and such that fi Sl i = fii and the total angular momentum is conserved. Since 
L*,i <C L\ for realistic parameters for a hierarchical triple system, the coupled tidal evolutions ensure that L\ ~ —£23- 

We also include orbital evolution due to induced eccentricity, magnetic braking, and gravitational radiation, as described 
in Appendix A. However, we find that the timescales associated with these processes are generally longer than the lifetime of 
a system such as HD 181068, so we do not discuss these effects in detail below. 



7.1 Results of Orbital Evolution 

Figure [8] shows an example of our evolutionary calculations. We plot the stellar radius, R\, orbital frequencies fii and fi23, 
and orbital decay time scales (ttide> t ccc , £ma g , and iGw) as a function of time. We begin our orbital calculations as Star 
1 is moving off the main sequence. The properties of the oscillation modes (calculated for helium burning red-giant model 
described in Section 6.1 1 are not appropriate for the initial model main sequence model, but become realistic as Star 1 moves 
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Figure 8. Top: Radius R\ as a function of time since leaving the main sequence for an M\ = 3.0Mq stellar model. The horizontal blue 
strip is the measured radius of HD 181068, within the measurement uncertainty. Middle: Angular orbital frequencies f!23 (solid black 
line) and 10f2i (solid red line) as a function of time, including tidal dissipation. The horizontal lines are the values of Q23 (dashed black) 
and 10f2i (dashed red) for the HD 181068 system. Bottom: tidal dissipation timescale t t jd e (solid black line), induced eccentricity orbital 
decay time scale t ccc (triple dot dashed green line), gravitational wave decay timescale £gw (dotted blue line), magnetic braking orbital 
decay time scale t ma g (dot dashed red line), and stellar evolution time scale tn (dashed purple line) as a function of time. 
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up the red giant branch. We begin our calculation with orbital frequencies of ^23,1 = 0.2D.23, o and = 1.15fli, , where 
^23,0 and Qi )0 are the observed orbital frequencies in HD 181068. 

Let us start by examining the top panel of Figure [8] The spike in radius at t — 1 x 10 7 yr corresponds to the largest radius 
obtained during the red giant phase, while the long flat period between 2 x 10 7 yr < t < 10 8 yr is the core helium burning 
phase. Figure [8] indicates that Star 1 in HD 181068 could be ascending the red giant branch, but is most likely a horizontal 
branch star. The measured temperature of Star 1 is also consistent with these possibilities. 

The middle panel of Figure [8] shows that the orbits evolve substantially due to tidal dissipation. The outer orbital 
frequency, £7i becomes slightly smaller as angular momentum is transferred from the inner to the outer orbit. The inner 
orbital frequency Q23 changes substantially, increasing by a factor of roughly five. The tidal dissipation begins as Star 1 moves 
up the red giant branch, increasing the value of i?i and hence Z/23- Tidal dissipation also becomes much stronger during this 
stage because the frequencies of envelope p-modes (which couple strongly with the tidal potential) become comparable to the 
tidal forcing frequencies. The tidal dissipation remains strong as Star 1 shrinks and descends the red giant branch because 
the system enters into a resonance locking configuration (see Section 7.2|. Once Star 1 settles onto the horizontal branch, the 
resonance locking ends and tidal dissipation is drastically reduced. 

The bottom panel of Figure [8] displays the relevant time scales of the evolution, including the stellar evolution time scale 
defined as 

*« = TiTT" (52) 

When Star 1 is on the main sequence, all the orbital time scales are longer than a Hubble time and can be ignored. As Star 
1 moves up the red giant branch, the tidal dissipation time scale ttide decreases by several orders of magnitude, as discussed 
above. When Star 1 is descending the red giant branch (at t « 1.5 x 10 7 yr), we see that tyde — in- This is a natural consequence 
of resonance locking, which we discuss in Section [7. 2| 

In contrast, the value of ttide is generally much larger during the core helium burning phase, due to the decreased values 
of _Ri and 0,23. While in the core helium burning phase, the value of ttide has sharp dips due to resonance crossings with mixed 
modes. There are also deeper, broader dips due to resonance crossings with envelope p-modes, such as the resonance crossing 
that occurs at t ~ 9 x 10 7 yr. However, these resonance crossing events are fairly brief and produce only small amounts of tidal 
dissipation. Thus, we are thus unlikely to observe any orbital decay in HD 181068, unless the system is in a resonance-locking 
state. 

The initial orbital frequencies were chosen such that the final values of Q23 and fii would be roughly equal to the 
observed orbital frequencies in HD 181068. This does not entail that our chosen values of ^23,1 and fii,* were the actual 
orbital frequencies of HD 181068 while Star 1 was on the main sequence. In reality, the initial orbital configuration of HD 
181068 may have contained inclined or eccentric orbits for which the tidal dissipation rates may be substantially different. 
Nonetheless, the fact that a small value of fl23,i is required to match observations of HD 181068 indicates that the system 
may have experienced substantial orbital decay due to three-body tidal dissipation. 

In the future, HD 181068 will likely enter a common envelope phase as Star 1 evolves up the asymptotic giant branch 
and envelops Stars 2 and 3. Additionally, it is possible that HD 181068 reached its compact configuration through a mass 
transfer phase while Star 1 was on the red giant branch. Our MESA-generated stellar models indicate that our 2>Mq stellar 
model would not have overflown its Roche lobe (assuming a\ was equal to its current value) while Star 1 was on the red giant 
branch. However, a 4Mq star probably would have overflowed its Roche lobe while on the red giant branch, and it may be 
possible that this occurred in HD 181068. However, the outcome of stable mass transfer or common envelope evolution in 
triple systems is highly uncertain, and it warrants further study. 



7.2 Resonance Locking 

The orbital evolution discussed above contains periods of resonance locking in which a tidally excited mode is held near 
resonance for long periods of time, causing large amounts of tidal dissipation. Resonance locking involving tidally driven 
stellar oscillation modes was first investigated by Witte & Savonije (1999) and recently proposed by Fuller & Lai (2012) to 
explain the observed tidally driven oscillations in KOI-54. Resonance locking can occur when the frequency uj a of a mode 
changes due to a stellar evolutionary process. As the mode approaches resonance with a tidal forcing frequency, enhanced 
tidal dissipation occurs. The dissipation feeds back on the process, changing the value of the forcing frequency Vf by spinning 
up the star or causing orbital decay. Under the right circumstances, the system maintains the nearly resonant configuration 
such that uJa — Vf, i- e -, the orbital evolution timescale is roughly equal to the stellar evolution time scale. 

The resonance locking shown in Section [7] is qualitatively different from the locking investigated by Witte & Savonije 
(1999) and Fuller & Lai (2012), although the general principles described above are still true. In the case of a system like HD 
181068, when Star 1 is descending the red giant branch, the frequencies of the stellar oscillation modes are increasing because 
the dynamical frequency of the star is increasing. When the value of ui a of an envelope mode is nearly resonant with a forcing 
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frequency Vf, enhanced tidal dissipation occurs, causing the orbit of Stars 2 and 3 about one another to decay. Hence, the 
values of 5^23 and Vf correspondingly increase, causing resonance locking. 

When the system is resonantly locked, uj a ~ uj. Assuming the star's oscillation frequencies change primarily due to the 
change in stellar radius, we have 

§£• (53) 
Furthermore, the tidal forcing frequency changes primarily due to the increasing value of ^23, so that 

^- 3 ^. (54) 

Vf 2fl 2 3 

Therefore, during resonance locking, 

W^IT' (55) 

Hi CL23 

and thus t t ido — tit. This explains the near equality of ttide and tit during the descent of Star 1 from the red giant branch in 
Figure [8] It also suggests that compact triples may endure a period of rapid tidal dissipation, caused by resonance locking, 
as the primary descends from the red giant branch toward the horizontal branch. Such tidal dissipation may lie in the near 
future for other observed compact hierarchical triples, such as KOI-126 (see Carter et al. 2011) and KOI-928 (see Steffen et 
al. 2011). 



8 DISCUSSION 

We have demonstrated that a new tidal dissipation mechanism exists for stars in compact hierarchical triple systems, even 
after the system has reached the quasi-equilibrium state of aligned and circular orbits with aligned and synchronous stellar 
spins. The three-body tidal forces produce forcing at frequencies of a = 2(^23 — 2fli), a = 2(^23 — fii), and a — 2Q2S 
(although more forcing frequencies will exist for non-coplanar and non-circular orbits) . If the primary star is a red giant with 
a large radius and thick convective envelope, the three-body tidal potential can couple strongly with the envelope p-modes of 
the star, exciting modes to large amplitudes. 

We compare our results to Kepler observations of HD 181068. The presence of oscillations in the lightcurve of HD 181068 
at frequencies a — 2(^23 — 2fii), a — 2(f223 — fii), and a — 2f223 indicates that the primary exhibits stellar oscillation modes 
excited by the three-body tidal forcing described in this paper. The large amplitude of the oscillations in HD 181068 are 
either due to large temperature variations produced by small amplitude, non-resonant modes, or they are due to geometric 
distortions produced by large amplitude, nearly resonant modes. In the latter case, the tidal dissipation rate is rapid, and the 
modes may be locked in resonance (see Section [7. 2\ . However, the amplitudes and phases of the modes are best explained by 
the non-resonant scenario. Furthermore, since the lifetime of the resonance locking phase is brief compared to the lifetime of 
the star on the horizontal branch, we find the non-resonant scenario to be much more likely. 

Because we have assumed adiabaticity when calculating stellar oscillation modes, there remains some uncertainty in the 
precise mode visibilities and damping rates. Fully non-adiabatic mode calculations can constrain the values of H a and 7^ in 
future studies. Furthermore, our calculations have been limited to the linear regime. As mentioned in Section[2] the three-body 
tidal forcing frequencies are nearly identical in the rotating frame of Star 1 if it is synchronized with the outer binary. This 
may allow for greatly enhanced non-linear mode coupling, affecting the mode amplitudes, damping rates, and visibilities. 

The three-body tidal effects can also produce substantial orbital evolution among the stellar components, the main effect 
of which is the decay of the inner orbit of the compact binary. Our orbital evolution calculations for a system resembling HD 
181068 reveal that, at most times, three-body tidal dissipation acts on long time scales and can be ignored. However, when the 
primary star is high on the red giant branch, three-body tidal dissipation may cause substantial orbital decay. Furthermore, 
stellar oscillation modes can become locked in resonance as the primary descends the red giant branch, resulting in greatly 
enhanced tidal dissipation. During resonance locking, tidal dissipation occurs on the same timescale as the stellar evolution, 
such that the orbital semi-major axis of Stars 2 and 3 decays as 023/023 — Ri/Ri- 

Future observations can detect and characterize tidally excited modes in compact triples. In non-eclipsing systems, the 
three-body nature of compact triples may not initially be detected from photometric and spectroscopic observations, especially 
if one star is much more luminous than its companions. The Kepler public red giant sample, consisting of more than 15,000 
stars (Hekker et al. 2011), provides a unique resource that may contain several hidden clones of the Trinity system, which 
could be used to test our theory in a broader parameter space. The signature of three-body tidal forcing is a triplet of evenly 
spaced modes with frequencies a — 2(^23 — 2fii), a = 2(^23 — fii), and a = 2Q.2Z, although the highest oscillation frequency 
may have a very small amplitude. The triplet could be mistaken for a rotationally split triplet of solar-like oscillations, but 
the tidally excited modes can be distinguished by their indefinitely long lifetimes (and thus narrow Fourier peaks). In low 
inclination systems, the a = 2(^23 — ^1) oscillation will produce the largest luminosity fluctuations, while the a = 2(^23 — 2fii) 
oscillation will be dominant in high-inclination systems like HD 181068. 



18 J. Fuller et al. 



ACKNOWLEDGMENTS 

We thank Dong Lai, Lars Bildsten, and Dan Tamayo for useful discussions. JF acknowledges the hospitality (Fall 2011) of 
the Kavli Institute for Theoretical Physics at UCSB (funded by the NSF through PHY 11-25915) where part of the work was 
carried out. This project has been supported by NSF grants AST-1008245 and AST-1211061, NASA grants NNX12AF85G and 
NNX10AP19G, the Hungarian OTKA Grants K76816, K83790 and MB08C 81013, and the "Lendiilet-2009" Young Researchers 
Program of the Hungarian Academy of Sciences. Part of the research leading to these results has received funding from 
the European Research Council under the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC 
grant agreement n° 227224 (PROSPERITY). AD gratefully acknowledges financial support from the Magyary Zoltan Public 
Foundation. AD was supported by the Hungarian Eotvos fellowship. AD has been supported by the Janos Bolyai Research 
Scholarship of the Hungarian Academy of Sciences. LLK wishes to thank support from the European Communitys Seventh 
Framework Programme (FP7/2007- 2013) under grant agreement number 269194. The Kepler Team and the Kepler Guest 
Observer Office are recognized for helping to make the mission and these data possible. 



REFERENCES 

Alexander, M.E., 1973, Astrophys. Space Sci., 23, 459 
Baker, N. & Kippenhahn, R., 1965, ApJ, 142, 868 
Beck, P., et al., 2012, Nature, 481, 55 

Bedding, T., et al., 2011, Nature, 471, 608 

Borkovits, T., et al., 2012, MNRAS, accepted | |arXiv:1210.1061| l 
Breger, M., et al., 1993, A&A, 271, 482 

Brown, T., Gilliland, R., Noyes, R., Ramsey, L., 1991, ApJ, 368, 599 
Burkart, J., Quataert, E., Arras, P., Weinberg, N. 2011, MNRAS, 421, 983 
Buta, R.J. & Smith, MA. 1979, ApJ, 232, 213 
Carter, J., et al. 2011, Science, 331, 562 
Chaplin, W., et al., 2011, Science, 332, 213 
Chaplin, W., et al. 2011, ApJL, 732, 5 

Correia, A., Laskar, J., Farago, F., Bou, G. 2011, CeMDA, 111, 105 

Derekas, A., et al., 2011, Science, 332, 216 

De Ridder, J., et al., 2009, Nature, 459, 398 

Dzicmbowski, W. 1971, Acta Astronomica, 21, 3 

Eggleton, P., Kisleva-Eggleton, L., 2001, ApJ, 562, 1012 

Fabrycky, D., Tremaine, S., 2007, ApJ, 669, 1298 

Ford, E., Kozinsky, B., Rasio, F., 2000, ApJ, 535, 385 

Fuller, J., Lai, D. 2012, MNRAS, 420, 3126 

Georgakarakos, N., 2002, MNRAS, 337, 559 

Gouttebroze, P., Toutain, T., 1994, A&A, 287, 535 

Hekker, S., et al., 2011, MNRAS, 414, 2594 

Houdek, G., Balmforth, N. J., Christensen-Dalsgaard, J., Gough, D. O., 1999, A&A, 351, 582 

Huber, D., et al., 2010, ApJ, 723, 1607 

Huber, D., et al., 2011, ApJ, 743, 143 

Hut, P. 1981, AA, 99, 126 

Kallinger, T., et al., 2010, A&A, 522, Al 

Kumar, P. & Quataert, E.J. 1998, ApJ 493, 412 

Lai, D. 1996, ApJ, 466, L35 

Lai, D. 1997, ApJ, 490, 847 

Lenz, P., & Breger, M., 2005, Comm. Asteroseis., 146, 53 
Lithwick, Y. & Wu, Y., 2012, ApJL, 756, 11 
Mathur, S., et al., 2011, ApJ, 741, 119 
Mazeh, T., Shaham, J., 1979, A&A, 77, 145 
Mosser, B., et al., 2012, |arXiv:1209.3336] 

de Medeiros, J.R., Da Rocha, C, Mayor, M., 1996, A&A 314, 499 

Paxton, B., Bildsten, L., Dotter, A., Herwig, F., Lasaffre, P. Timmes, F. 2011, ApJS, 192, 3 
Robinson, E., Kepler, S., Nather, R. 1982, ApJ, 259, 219 
Samadi, et al., 2007, A&A, 463, 297 

Savitzky, A., Golay, M., 1964, Analytical Chemistry, 36, 1627 

Schenck, A. K., Arras, P., Flanagan, E. E., Teukolsky, S. A., Wasserman, I. 2001, Phys Rev D, 65, 024001 
Steffen, J., et al. 2011, MNRAS, 417L, 31 
Thompson, S., et al. 2012, ApJ, 753, 86 

Unno, W., Osaki, Y., Ando, H., Saio, H., Shibahashi, H. 1989, Nonradial Oscillations of Stars (University of Tokyo Press) 

Welsh, W. F., et al. 2011, ApJS, 197, 4 

Wittc, M. G. & Savonije, G.J. 1999, ApJ, 350, 129 



Tides in Triple Systems 19 



APPENDIX A: NON-TIDAL ORBITAL EVOLUTION 

We wish to compare the effects of three-body tides to other physical mechanisms that may produce orbital evolution in 
triple-star systems. One such mechanism is tidal orbital decay due to the induced eccentricity of the inner binary (Stars 2 
and 3) due to the perturbing gravitational influence of Star 1. This effect has been studied in detail by, e.g., Ford et al. 2000. 
In this mechanism, the eccentricity of the orbit of Stars 2 and 3 about one another oscillates around a small, non-zero value 
due to the presence of Star 1. Tidal interactions between Stars 2 and 3 act to re-circularize the orbit, dissipating energy and 
causing the orbit of Stars 2 and 3 about one another to decay. The orbital decay time scale due to dissipation in Star 2 is 
(Lithwick & Wu 2012) 

023 = 1 Q O23-B2 1 (a2 3 \ 8 1 



ccc ~ h^i ~ 18^~h~GM^ q(l + q) \Jh) (P)' 1 ' 

where Q is the tidal quality factor, k 2 is the constant of apsidal motion, q = M2/M3, and (e 2 ) is the mean square eccentricity 
of the orbit of Stars 2 and 3 about one another. We have divided by an extra factor of two to account for tidal dissipation in 
both stars. According to Georgakarakos (2002), the mean square eccentricity for a nearly equal mass binary (M 2 ~ M3) is 

However, for unequal mass inner binaries (M 2 7^ M3) or systems near an orbital resonance, the mean square eccentricity may 
be significantly larger. In our calculations, we use k 2 — 0.02 (appropriate for low mass dwarf stars), and Q = 10 5 . 

Another effect that can cause orbital decay in a compact triple system is magnetic braking. Magnetic braking is likely 
to act on short time scales in systems like HD 181068 because Stars 2 and 3 are rapidly rotating, low mass stars with 
convective envelopes that likely have strong magnetic dynamos. We estimate the torque due to magnetic braking according 
to the prescription of Krishnamurthi et al. 1997, yielding a magnetic braking timescale of 

a 23 ^iah f m 2 \ i/2 f R 2 y 1/2 

tmag a 2:i AK MB ni c \M Q ) \RcJ ' 

where f2 SjC is a critical rotation rate of Q. c ~ lOfi Sj0 , and Kmb is a calibrated constant of Kmb ~ 2.6 x 10 47 g s cm 2 . Once 
again, we have divided by an extra factor of two to account for magnetic braking produced by Stars 2 and 3, assuming they 
are nearly identical. Star 1 may produce additional magnetic braking, but because a\ ^> 023 in hierarchical triples, the time 
scale of tidal orbital decay associated with this process is long, and can be ignored. 

Finally, the orbit of Stars 2 and 3 may decay due to the emission of gravitational waves. The orbital decay time scale due 
to gravitational waves is 

= 023 = J3C 5 O23 (A4) 

GW a 23 64G 3 M 2 M 3 (M 2 + M 3 ) ' 1 ' 



(A3) 



